Efficient implicit algorithm for the equations of 2-D viscous compressible flow: application to shock-boundary layer interaction

A fully implicit algorithm has been developed to time integrate the equations of 2-D compressible viscous flow. The algorithm was constructed so as to optimize computational efficiency. The time-consuming block matrix inversions usually associated with implicit algorithms have been reduced to the trivial non-iterative inversion of four sets of scalar bidiagonal matrices. Thus, the algorithm requires virtually no more computer storage than an explicit algorithm. The efficient structure of the implicit algorithm is reflected in comparative timings which slow that it requires only a factor of two more computer time per point per time step than a typical explicit algorithm. Therefore, the algorithm allows more economical solution of given flows than existing explicit methods and also allows more difficult problems to be attempted using available computer resources. Application of the algorithm to the problem of shock-boundary layer interaction produces results consistent with both experimental measurements and other calculations.     

A Multibody Loop Constraints Approach for Modelling Cam/Follower Devices

This paper proposes a formulation for modelling mechanisms witha cam/follower type of contact using a multibody approach in relativecoordinates. The proposed approach is inspired from the wheel/railcontact model developed by Fisette and Samin but, in the present case,possible intermittent contact between the cam and the follower isconsidered, for generality purposes.Loop kinematicconstraints are introduced to satisfy tangent and punctual contact aslong as the bodies lean against each other. The effective presence (ornot) of the contact is governed by the sign of the normal constraintforce which can be computed thanks to the Lagrange multiplierstechnique.The above-mentioned option to kinematicallyconstraint the bodies in their contact phase unavoidably leads to ashift from one model to another when a contact disappears (or,conversely, reappears). Indeed, this increases (or decreases) the numberof degrees of freedom of the current system. The control of the variablepartitioning is thus absolutely necessary and is all the more complex sothat practical applications can contain several pairs of bodies inintermittent contact.Illustrative examples are proposed: acomparison with another multibody formalism, an experimental validationand the modelling of universal wheels of mobile robots which representsa quite original application of the proposed formulation.     

卷积型加权残值法求解薄壳的动力学问题

板壳等弹性体,受到动荷载作用时,其动力分析问题是比较重要且难于解决的．利用卷积型加权残值法,推导出Gurtin变分原理,并应用卷积型加权残值法计算了薄壳的动力学问题,为计算薄壳的结构动力学问题提供了一种有效的方法．     

Nonlinear Responses of a Magnetoelastic Beam in a Step-Pulsed Magnetic Field

A theoretical study is carried out on the dynamics of a magnetoelastic beam being in a step-pulsed magnetic field. For this aim, the magnetic potential and elastic energies are determined for the beam and partial differential equations are established according to Hamilton's principle. It is proven that the magnetoelastic beam can give a variety of complex behavior in the case of step-pulsed field excitations. An intermediate regime of two-well chaos is observed. Theoretical findings were found to be in a good agreement with the experimental results for the specific system parameters. On leave from Institute of Physics, University of Bayreuth, 65440 Bayreuth, Germany An erratum to this article is available at .     

Almost Automorphic and Almost Periodic Dynamics for Quasimonotone Non-Autonomous Functional Differential Equations

The occurrence of almost automorphic dynamics for monotone non-autonomous recurrent finite-delay functional differential equations is analyzed. Topological methods are used to ensure its presence in the case of existence of semicontinuous semi-equilibria. When these semi-equilibria are continuous and strong, the presence of almost automorphic extensions is persistent under small perturbations. The above method provides a minimal set isomorphic to the base in the case of a convex semiflow. Some examples show the applicability of these results.     

Stability and Oscillations in a Time-Delayed Vehicle System with Driver Control

In this paper, linear stability and chaotic motion of a time-delayednonlinear vehicle system are studied. The stability is determined bycomputing the spectrum associated with a system of linear retardedfunctional differential equations, which reveals that a loss ofstability occurs following a Hopf bifurcation. Beyond the critical valuefor linear stability, the system exhibits limit cycle motions.Subharmonic, quasi-periodic and chaotic motions are observed for asystem excited by a periodic disturbance.     

Non-linear normal modes,invariance, and modal dynamics approximations of non-linear systems

Non-linear systems are here tackled in a manner directly inherited from linear ones, that is, by using proper normal modes of motion. These are defined in terms of invariant manifolds in the system's phase space, on which the uncoupled system dynamics can be studied. Two different methodologies which were previously developed to derive the non-linear normal modes of continuous systems — one based on a purely continuous approach, and one based on a discretized approach to which the theory developed for discrete systems can be applied-are simultaneously applied to the same study case-an Euler-Bernoulli beam constrained by a non-linear spring-and compared as regards accuracy and reliability. Numerical simulations of pure non-linear modal motions are performed using these approaches, and compared to simulations of equations obtained by a classical projection onto the linear modes. The invariance properties of the non-linear normal modes are demonstrated, and it is also found that, for a pure non-linear modal motion, the invariant manifold approach achieves the same accuracy as that obtained using several linear normal modes, but with significantly reduced computational cost. This is mainly due to the possibility of obtaining high-order accuracy in the dynamics by solving only one non-linear ordinary differential equation.     

Direct modeling of flow of FENE fluids

Direct simulations of macromolecular fluids are carried out for flows between parallel plates and in expanding and contracting channels. The macromolecules are modeled as FENE dumbbells with soft disks or Lennard-Jones dumbbell-dumbbell interactions. The results are presented in terms of profiles and contour plots of velocity, pressure, temperature, density, and flow fields. In addition the data for potential energy, shear stress, and the normal components of the stress tensor are collected. In general, an excellent agreement is found between the simulated profiles and the well-known flow structures, such as flow separation and formation of viscous eddies, indicating that micro-hydrodynamics is a viable tool in linking macroscopic phenomena with the underlying physical mechanisms. The simulations are performed in the Newtonian regime, for medium-size systems comprising up to 3888 dumbbells. This number is sufficiently large to control boundary and particle number effects. The flow is induced by gravity. The traditional stochastic (thermal) and periodic boundary conditions are employed. Also, diffusive boundary conditions, which could include a stagnant fluid layer and repulsive potential walls, are developed. The scaling problems, which are related to the application of a large external force in a microscopic system (of the size of the order 100 Å), result in extreme pressure and temperature gradients. In addition, the viscosity and thermal conductivity coefficients obtained from velocity and temperature profiles of the channel flow are presented. These results are confirmed independently from modeling of Couette flow by the SLLOD equations of motion and from the Evans algorithm for thermal conductivity.     

Space-time perturbation modes for non-linear dynamic analysis of beams

In this work an attempt is made at bridging the powerful perturbation methods of analytical dynamics to the versatile finite element techniques which can readily handle arbitrarily complex structures. The proposed analysis methodology has two distinguishing features. First, a space-time finite element formulation is used, and hence the concept of modes is here naturally extended to that of space-time modes, where the time dependency is implied in the assumed modes. As a result, the partial differential equations of motion are directly reduced to purely algebraic non-linear simultaneous equations. Second, perturbation modes, rather than the usual vibration mode shapes are used and shown to be an appropriate basis for non-linear dynamic analysis. These modes bring information about the non-linearities of the system through the higher order derivatives of the strain and kinetic energies. The procedure is illustrated on non-linear beam problems and the results are compared with those of a full finite element model, i.e., when all the degrees of freedom are considered, as well as with analytical results, when available.     

Numerical Analysis of the Influence of Thermal Stresses on the Nonlinear Thermomolecular Pressure Difference Effect

On the basis of a numerical analysis of the non-Navier-Stokes gas-dynamic equations for slow non-isothermal gas flows, the nonlinear thermomolecular pressure difference effect due to a large temperature gradient along the lateral surface of a capillary is investigated. It is shown that the magnitude of the effect is substantially different from the values calculated using the Navier-Stokes equations. For two models of molecular interaction (Maxwell molecules and hard spheres), the possibility of a quasi-one-dimensional interpretation of the effect for experimental estimation purposes is demonstrated. The solutions of the relaxation kinetic equation for flow in a plane capillary at small Knudsen numbers and the gas-dynamic equations for slow non-isothermal flows are compared and the range of their applicability is estimated.